Inhibiting Three-Body Recombination in Atomic Bose-Einstein Condensates 
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We discuss the possibility of inhibiting three-body recombination in atomic Bose-Einstein con- 
densates via the application of resonant 2ix laser pulses. These pulses result in the periodic change 
in the phase of the molecular state by 7r, which leads to destructive interference between the decay 
amplitudes following successive pulses. We show that the decay rate can be reduced by several 
orders of magnitude under realistic conditions. 
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A fundamental limit to the lifetime of trapped atomic 
Bose-Einstein condensates (BECs) is three-body recom- 
bination, in which two atoms form a molecular dimer 
and the third atom carries away the excess energy and 
momentum Q, 0, EJ El ■ Unlike other loss mechanisms 
such as spin relaxation, collisions with background ther- 
mal atoms, or spontaneous light scattering, three-body 
recombination is an intrinsic loss mechanism that cannot 
be eliminated by simply engineering a better trapping 
environment. 

The decay rate due to three-body recombination is pro- 
portional to the square of the atomic density, 



in 



where n is the atomic density, K3 ~ ha 4 /m, and a is 
the elastic atomic s-wave scattering length for atoms of 

mass m S IE Q, B H E3- ( We note that E q- © 
contains an additional factor of 1/6 if all atoms occupy 
the same quantum state Q). At high densities and large 
scattering lengths, three-body recombination represents 
the primary limit on the lifetime of condensates. 

The strong dependence of three-body recombination 
on a and the density imposes severe restrictions on ex- 
periments done in the vicinity of a Feshbach resonance 
HHH El, where a diverges, as well as for tightly confined 
samples. The latter situation is relevant in particular 
to the new field of integrated atom optics 0] . Thus it 
is important to pursue possible techniques for inhibiting 
three-body recombination in condensates. 

The control of decay and decoherence mechanisms in 
quantu m sy stems is also important in other areas of 
physics |l2t 1 1 3| . Much of this work is motivated by the 
requirements for quantum computers as well as an in- 
terest in fundamental issues in quantum mechanics, such 
as the Quantum Zeno effect [l^. An important result 



in this context is the demonstration that time-dependent 
external fields that modulate either the energies of the 
system under consideration or its coupling to the exter- 
nal reservoir into which it decays can reduce its decay 
rate. A key requirement of such schemes is that the ex- 
ternal field be modulated on a time scale shorter than 
the reservoir correlation time. 

For positive scattering lengths much larger than the 
range of the interatomic potential, three-body recom- 
bination is dominated by the formation of molecules in 
a weakly bound state with binding energy £{, = Ti/ma 2 
2, 3]. In this paper we discuss the use of a regular se- 
quence of resonant 2ir laser pulses applied to that state 
to inhibit three-body decay. The effect of each pulse is 
to change its phase by ir. This leads to a destructive 
interference in the probability amplitudes between suc- 
cessive pulses for the formation of a molecule. A similar 
idea was explored by Agarwal et al. for inhibiting spon- 
taneous emission from a two-level atom |l3j . However, it 
is virtually impossible to inhibit spontaneous emission in 
free space because of the nearly instantaneous correlation 
time of the vacuum. The situation is more favorable in 
the present case: The role of the reservoir is now played 
by the molecular states formed by three-body recombi- 
nation, and its correlation time is roughly on the order 



of the inverse of the binding energy, e b 
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for 



a ~ lOOOao and to ~ 50 a.m.u. Therefore a sequence 
of pulses separated by an interval of T <C e^ 1 should be 
effective at decreasing the rate of molecule formation. 

We proceed by deriving a coarse-grained master equa- 
tion for the evolution of a BEC in the presence of a se- 
quence of impulsive laser pulses. The pulses are assumed 
to be off-resonant with respect to the condensate atoms, 
but resonant with respect to the molecules formed by 
three-body recombination. The total Hamiltonian for the 
atom-molecule system is H = H a + H m + H3 where 
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Here 7J a and i? m are the atomic and molecular Hamil- 
tonians, respectively, and H 3 accounts for three-body re- 
combination yj. The field operators ^>(r) and </> g (r) de- 
scribe the annihilation of atoms and of molecules with 
binding energy e&, respectively. In H a , g = A-Kfi 2 a/m 
while V(r) is an external trapping potential and p is the 
chemical potential. 

We assume that the laser field, with Rabi frequency 
fi m (i) and frequency wl, is resonant with a single vi- 
brational state in the molecular potential of the electron- 
ically excited molecule. This state is denoted by the anni- 
hilation operator <p e (r) . Its binding energy relative to the 
electronic energy of the corresponding free atoms is e e , so 



that the resonance condition is 5 = (uj e 



,) + e b -LU L 



0. We further assume that the laser is off- resonant with 
respect to the excited state transition for the correspond- 
ing free atoms, A a = u eg -wj, > fl™ / y/2, 7„, , where 
ttm/VZ is the atomic Rabi frequency [15j and j a the 
linewidth. For 6 = this gives A a = e e — e^. Since e& 
corresponds to a weakly bound state near the dissocia- 
tion limit, e e would have to correspond to a low lying vi- 
brational state in the molecular potential for the excited 
atoms with a binding energy of the order ~ 1 — 10GHz. 

For the problem at hand, f2 m (t) = Clof(t) where f(t) is 
a train of square pulses of unit amplitude with duration 
t p and separation T such that £Iqt p = 2-7T. If r p is short 
compared to the characteristic time of the center-of-mass 
molecular dynamics, the molecules undergo a complete 
Rabi oscillation for each pulse, leaving their excited state 
unchanged while their ground state acquires a ir phase 
shift. We can then eliminate the excited molecular state 
from H m and make for the ground state field operators 
the substitution 



ff (r, *)-»(-!) 



-1) [t+Tp\ 



where [...] denotes the integer part of the term in brack- 
ets. 

The Hamiltonian H3 describes the formation of dimers 
from three colliding atoms in which the third atom car- 
ries away the excess kinetic energy and momentum re- 
leased by the molecule formation. The binding energy of 
the molecule is converted into molecular and atomic ki- 
netic energies, t\, = 3fiK 2 /Am where +KK and — fiK are 
the atomic and molecular momenta. For fiK 2 /4m 3> p, 
we can decompose the atomic field operator as tp(r) = 
ipT(r) + V'F(r) where %1>t represent the trapped conden- 
sate atoms while i\}p are atoms with kinetic energy 2ef,/3. 
By making what amounts to the rotating wave approxi- 
mation and keeping only resonant terms, H3 reduces to 



H 3 = Kk I d 3 xft^pipT + h- c - 

We note that the molecules and recoiling atoms es- 
cape from the trap for e&/3 > Vo where Vb is the trap 
depth. For large scattering lengths with e& ~ 10 s , 
atoms and molecules are lost for Vb/fcs < WpK. Un- 
der these conditions, it is convenient to adopt a plane 
wave basis for the atoms and molecules that are lost 
from the trap, (f> g (r,t) = (1/y/V) J2k &k(t) exp[ik ■ r] 
and ^ir(r,t) = (l/y/V) J2 k Ck(t) exp[ik • r], where V is 
a quantization volume. For the condensate we use a 
zero-temperature single- mode approximation, ipr(T, t) = 
uo(r)b(t) and Uq(t) is a Hartree wave function for the 
condensate ground state with eigenvalue p. After trans- 
forming to the interaction representation, the interaction 
Hamiltonian becomes 

H 3 (t) = ?iC/(ki+k 2 )(-l)[*/ T +^e^ 12t aj ii cj i2 6 3 +/i.c 



ki,k2 



(2) 



where Suji2 = /ikf/(4m) + /ik|/(2m) — th and U(k) = 
k J d z xe~ lk ' T ul(r)/V . 

Starting from Eq. J5J) we can derive a master equation 
for the density operator of the condensate atoms, p(t). 
In this approach the molecules and free atoms are inter- 
preted as a reservoir coupled to the condensate by H 3 (t) . 
It is important to note that H 3 (t) is not a continuous 
function of time since the pulses result in a discontinuous 
sign change after each interval T + t p . Hence one must 
proceed with caution when deriving the master equation 
since differentiation and integration are no longer inverses 
of each other. We use the form of the master equation 
derived e.g. in Ref. for the coarse-grained derivative 
p(t) = (p(t + r) — p(t))/r of the system density operator, 
r being a time long compared to the correlation time of 
the reservoir t c but short compared to times over which 
the condensate evolves, r c <C t <C I/711 2 . Taking the 
reservoir density operator to be in the vacuum state and 
assuming that r = NT where N is the number of pulses 
and r p <C T, we obtain 



vlH 3 lP (t) 

(S t3 S 3 p(t) + p{tp 3 b 3 - 2b 3 p{0 3 ^j . (3) 



Here v = is(T, N) is a "Lamb shift" of the condensate 
atoms. We do not reproduce its lengthy expression, con- 
centrating instead on the decay rate of the condensate, 
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7 = j(T, N) due to molecule formation, 

j(T,N) = |C^(ki +k 2 )| 2 tan 2 (^ 12 T/2) 



ki,k 2 



sin 2 (5lo 12 NT/2 + Nir/2) 
NT5uo 2 l2 /A 



(4) 



The tan 2 (&j 12 T/2) term in j(T, N) describes the ef- 
fect of the pulses on the decay of the condensate. We 
note that 7(T, N)NT agrees with the transition prob- 
ability calculated directly from H 3 (t) using first-order 
time-dependent perturbation theory |13| . 

For a Markovian reservoir, r c — > 0, one can let r = 
NT —>■ oo 0| . At the same time we note that for an even 
number of pulses no net phase (modulo 2tt) is acquired 
by the molecules. Under these conditions the decay rate 
is 



lim 7 = 2ir 



E 

ki,k 2 



\UQli + k 2 )| 2 tan 2 (6u 12 T/2) 6{6w 12 ). 



(5) 

For tan 2 (5lo 12 T/2) = 1, one then recovers the standard 
result for the decay rate in the Markov limit, as expected. 

It is clear from Eq. (|SJ that lim^oo 7 = 0. This is due 
to the fact that for 5uoi 2 — 0, the phases of the molecules 
and free atoms do not change in the interval between 
pulses, j(T + t p ) < t < (J + 1)T, where j is an integer. 
However at the end of each interval, r = j(T + t p ), the 
phase of the molecules changes by n. As a result, there is 
complete destructive interference between the transition 
amplitudes for neighboring intervals. 

For a finite reservoir correlation time r c , however, it is 
no longer possible to (formally) let r = NT — > 00 and 
the interferences cease to be fully destructive, resulting 
in a non-zero condensate decay rate. To reduce 7 be- 
low its unperturbed value requires tan 2 (5wi 2 T/2) <C 1 
but since r c is approximately given by the reciprocal of 
the reservoir bandwidth, e c , one must choose T -C 7rr c , 
that is, the laser pulses must be separated by an interval 
shorter than the correlation time of the molecular field. 
The bandwidth of the molecular reservoir is determined 
by the short-wavelength cutoff in which is on the 
order of the range of the interatomic potential, Rq < a 
[H [13 ■ This gives a bandwidth of e c ~ hjmR\ > tb- 

The de Broglie wavelengths of the molecules and en- 
ergetically free atoms are of the order a. For a trapped 
condensate of size £, one typically has I ^S> a. In this case, 
we may treat the gas as being locally homogeneous with 
respect to the fast moving molecules and atoms formed 
by recombination. For a uniform BEC, 7 then becomes 



7 = l0 F(Te b , N) 



(6) 



where 70 = (l/27r)(K/U) 2 (4m/3?i) 3/2 ^/e6 is the decay 
rate in the absence of any pulses in the t — > 00 Markov 
limit. Note that 7 is related to K 3 by K 3 = 367 0. The 



FIG. 1: Normalized decay rate for N = 20 and a — 1.5 
with pulses (solid line) and without pulses (dotted line) as a 
function of Teh 



dimensionless factor reflecting the effect of the 2it pulses 
is 



F(t,N) 



-taa 2 (xt/2) sin 2 ((a;t + tt)7V/2) 



27r7Vt(x/2) 2 



(7) 

Here a = (e c — eb)/tb is a dimensionless high energy cut- 
off. By setting tan 2 (xi/2) = 1 in Eq. (J7J we obtain 
the finite time decay rate in the absence of any pulses, 
which we denote as F?jp(t, N). F^p rapidly approaches 
its limiting value of 1 for t ~ 1 for all values of N and a 
considered below. This is illustrated in Fig. 1 for N = 20 
and a = 1.5 

The effective number of pulses N — [t/T] that con- 
tribute to the reduction of 7 is determined by the time 
over which the phases of the molecular states can evolve 
coherently. Any event that leads either to the decay of 
the molecular state or to a randomization of its phase 
will negate the accumulated effect of the pulses. Since 
the molecules formed by three-body recombination are 
in a very weakly bound state, they can decay to more 
deeply bound vibrational states via inelastic collisions 
with atoms [Isl flfll Efl l2ll | . The decay rate for the molec- 
ular state is then given by nn — f3nva ~ Anfinha/ySm. 
Here, nva is the elastic collision rate between atoms and 
molecules with velocity v = h/y/3ma and we assume that 
the cross section can be approximated by the atom-atom 
elastic cross section, a w 47ra 2 . f3 is then the ratio of the 
inelastic to elastic cross sections for transitions to deeply 
bound vibrational states of the molecules. 

We can therefore take r w l/ttn in order to determine 
the effective number of pulses. In this case the condition 
on the time scales involved in the derivation of the master 
equation, I/777 2 > I/7071 2 ^> r 3> r c , can be reexpressed 
as na 3 < 4nf3/^/3, V3/(47r/3)(a/i? ) 2 where na 3 is the 
dilute gas parameter and is typically <C 1. Refs. 0, 
Hil I20I l2l| give empirical values for k of 10~ 9 — 10~ n 
cm 3 /s in the vicinity of a Feshbach resonance. Using the 
values a ~ lOOOao and m ~ 50a. m.u. gives a range for 
P of 0.01 — 1. With these values of (3 one can reduce the 
above condition on the time scales for a > Rq to simply 
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FIG. 2: Normalized decay rate as a function of the dimen- 
sionless pulse period, Tet for N = 10 (solid line), N — 11 
(dashed line), N = 19 (dotted lined), and N = 20 (dashed 
dot) for a cutoff of a — 1.5. 
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neglected without loss of generality. Of greater concern 
are the atom losses due to Rayleigh scattering of laser 
photons. They should be less than the losses due to 
three-body recombination, otherwise we have simply re- 
placed one loss mechanism with another. For t 3> T + t p , 
the time-averaged Rayleigh scattering loss rate is f = 
( la W/ml)(r p /T) « (0.785/T)(|ft |/|A a |)( 7 „/|A a |) 
where we have again assumed that T ^> t p = l/27rfio- If 
we take j a ~ 10 7 s _1 , a detuning of A Q ~ 10 10 s _1 , and 
10r p = T ~ 10~ 6 s, we obtain a lifetime of f _1 ~ 10s. 
Note, however, that f is a single-particle loss rate and is 
unaffected by changes in the density n and a. 

In conclusion, we have discussed a technique for in- 
hibiting three-body recombination in BEC via a sequence 
of resonant 2ir laser pulses. We have shown that the 
three-body decay rate can be reduced to only a few per- 
cent of its value in the absence of the pulses. This method 
should be useful for extending the lifetime of condensates 
in the high density regime and near Feshbach resonances. 

This work is supported in part by the US Office of 
Naval Research, by the NSF, by the US Army Research 
Office, by NASA, and by the Joint Services Optics Pro- 
gram. 



FIG. 3: Normalized decay rate as a function of the dimen- 
sionless pulse period, Teb for N = 20 and a cutoff of a — 0.9 
(solid line), a — 1.3 (dashed line), a — 1.7 (dashed dot), and 
a = 2.1 (dotted lined). 

na 3 <C 0.1 so that the condensate must simply be dilute. 

These considerations imply that for a condensate den- 
sity of n = 10 15 /cm 3 and k « 10 _10 cm 3 /s, one can have 
on the order of 10 pulses within the lifetime of a molecule 
for a pulse period of T = 10~ 6 s. Fig. 2 shows a plot of 
7/70 as a function of T for N = 10, 11, 19, 20. It is clear 
that the more pulses that can be applied, the greater the 
reduction in the decay rate. In general for NTeb <§C 1 and 
finite, an even number of pulses, 2j, produces a lower de- 
cay rate than 2j ± 1 while for larger Teb the decay rate 
can be lower for an odd number of pulses. However, the 
difference between an even and odd number of pulses 
decreases as N gets larger. We note that F(t, N) is sen- 
sitive to the precise value of the cutoff 0] , as illustrated 
in Fig. 3. This is a common feature in effective low- 
energy field theories 0. Despite the dependence on a, 
one has F(Te b ,N) < 1 provided Te b ,Te b a < tt. Fig. 3 
shows that when Tt^a ~ it, the decay rate starts to in- 
crease rapidly. From these results, we conclude that the 
pulse train is able to reduce the decay rate to only a few 
percent of its "bare" value. 

Finally we discuss the effect of the laser pulses on the 
atoms. The atoms experience a periodic AC Stark shift, 
Ti\n a (t)\ 2 /4A a . It leads to a renormalization of the bind- 
ing energy relative to the dissociation limit and can be 
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